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Abstract 

We study the behavior of single linear polyelectrolytes condensed by trivalent salt under the 
action of electric fields through computer simulations. The chain is unfolded when the strength of 
the electric field is stronger than a critical value. This critical electric field follows a scaling law 
against chain length and the exponent of the scaling law is —0.77(1), smaller than the theoretical 
prediction, — 3z^/2 [Netz, Phys. Rev. Lett. 90 (2003) 128104], and the one obtained by simulations 
in tetravalent salt solutions, -0.453(3) [Hsiao and Wu, J. Phys. Chem. B 112 (2008) 13179]. It 
demonstrates that the scaling exponent depends sensitively on the salt valence. Hence, it is easier 
to unfold chains condensed by multivalent salt of smaller valence. Moreover, the absolute value of 
chain electrophoretic mobility increases drastically when the chain is unfolded in an electric field. 
The dependence of the mobility on electric field and chain length provides a plausible way to impart 
chain-length dependence in free-solution electrophoresis via chain unfolding transition induced by 
electric fields. Finally, we show that, in addition to an elongated structure, a condensed chain can 
be unfolded into an U-shaped structure. The formation of this structure in our study is purely a 
result of the electric polarization, but not of the elasto-hydrodynamics dominated in sedimentation 
of polymers. 
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I. INTRODUCTION 



To well understand the properties of charged macromolecules in electric fields, including 
the conformation and the mobility, is very important in many domains of researches such 
as in polymer science, in biophysics, and in microfiuidics, for the reason of a large variety of 
applications [1] . Applying electric fields stays at the center of the techniques to manipulate 
charged macromolecules. It can be also used as a tool to separate molecules by sizes. 
However, for the latter case, experiments are usually performed in a sieving matrix such 
as in gel, instead of in a solution This is because the free draining effect in an 

electrolyte solution produces an electrophoretic mobility independent of the chain length 
of macromolecules under a typical electrophoretic condition ^. Nevertheless, researchers 
continue to devote their efforts in finding ways for size separation in free solutions for the 
reason of its high throughput and applications in microfiuidics. 

In 2003, Netz proposed a new strategy to achieve this goal by unfolding condensed poly- 
electrolytes (PEs) in electric fields jsl, . He predicted that the chain mobility increases when 
a condensed PE chain unfolds in an electric field, and the critical electric field to unfold a 
chain, E*, depends on the chain length N, following the scaling law E* ~ A^^3i//2 ^j^g^-g jy 
is the chain swelling exponent. Therefore, longer chains will be unfolded and separated out 
earlier when the applying electric field slowly increases. His idea has been recently verified 
by simulations |3] in which PE chains were condensed into globules by tetravalent salts and 
then stretched in electric fields. A more general form of the scaling law has been proposed 
in the study, reading as E* ~ V^^^"^ where V is the ellipsoidal volume calculated from the 
three eigenvalues of the chain gyration tensor. According to the scaling law obtained by the 
simulations, an electric field of 2kV/cm should be applied to unfold collapsed PEs of chain 
length of order 10^. This electric field is relatively strong. 

For practical reason, we wish E* to be as small as possible. One way to reduce E* is to 
make less compact the condensed chain structure. This aim can be achieved, for example, by 
increasing temperature, by performing experiments in high-dielectric solutions, or by using 
weak condensing agent to collapse PEs. In this paper, we choose the last method, using 
trivalent salt as the condensing agent, and study the static and dynamic properties of chains 
and the unfolding electric field. Since the electrostatic interaction with trivalent salt is 25% 
weaker than that with tetravalent salt, some unexpected situations may take place. A key 
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question is to know if E* still follow the same scaling law of the strong condensation as 
shown in Netz's js, 6| and in Hsiao and Wu's study Q]- The rest of this paper is organized 
as follows. In Section II, we describe our model and simulation setup. In Section III, we 
present our results. The discussed topics include the degree of unfolding, the critical electric 
field to unfold a chain, the electrophoretic mobility, the distribution of condensed trivalent 
counterions on a chain, and the chain conformation after unfolding. We give our conclusions 
in Section IV. 

II. MODEL AND SIMULATION SETUP 

Our simulation system contains a single polyelectrolyte and trivalent salt, placed in a 
rectangular box with periodic boundary condition. The polyelectrolyte dissociates into a 
polyion chain and many counterions. The polyion is modeled by a bead-spring chain, con- 
sisting of beads; each bead carries a — e charge where e is the elementary charge unit. 
The counterions are modeled by spheres; each carries +e charge. The trivalent salt disso- 
ciates into trivalent cations (counterions) and monovalent anions (colons); these ions are 
also modeled by charged spheres. Solvent is treated as a uniform dielectric medium with 
dielectric constant equal to e^. Three kinds of interaction are considered: the excluded vol- 
ume interaction, the Coulomb interaction, and the bond connectivity. The excluded volume 
interaction is modeled by a purely repulsive Lennard- Jones potential 

, , / 4eLj [{cr/rf - {a/rf] + e,j for r < 2^^ 
Uex[r) = < (1) 
I for r > 2i/V 

where r is the distance between two particles, Slj is the interaction strength, and a denotes 
the diameter of a particle. We assumed that all the beads and spheres have identical Su and 
cr. We set Elj = ksT /1.2 where /c^ is the Boltzmann constant and T is the temperature. 
The Coulomb interaction is 

f/w(r) = (2) 

where Zj and Zj are the valences of the two charges and = / {Anere^kBT) is the 
Bjerrum length, at which two unit charges have the Coulomb interaction tantamount to 
the thermal energy /c^T. We set to be 3cr to simulate highly charged PEs, such as 
polystyrene sulfonate. Ucoui was calculated by PPPM Ewald method. Two adjacent beads 
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(monomers) on the chain are connected by the bond connectivity, modeled by a finitely 
extensible nonlinear elastic potential 

Ubondib) = -^hb'iax In - (3) 

\ max / 

where b is the bond length, b^ax is the maximum bond extention, and kh is the spring 
constant. We set bmax = 2(T and kf, = S.SSSS/csT/cr^. The average bond length under this 
setup is about l.lcx. An external uniform electric field E is applied, toward x direction. The 
equation of motion of a particle is described by the Langevin equation: 

dU 

miTi = -— - mi%fi + ZicEx + fji{t) (4) 

where is the mass of the particle rj is its position vector, mj7j is the friction coefficient, 
and 7]i simulates the random collision by solvent molecules. 'f]i{t) has zero mean over time 
and satisfies the fiuctuation-dissipation theorem: 

{Ut) ■ Vjit')) = 6kBTmaiSij5{t - t') (5) 

where 6ij and 6{t — t') are the Kronecker and the Dirac delta function, respectively. The tem- 
perature control is incorporated according to this theorem. We assumed that the particles 
have the same mass m and damping constant 7. We set 7 = Ir^^ where r = a^mj (ksT) 
is the time unit. We know that the dynamics of polymers in dilute solutions is described 
by Zimm model jsl. However, when an electric field is applied in a typical electrophoretic 
condition, the hydrodynamic interaction is largely canceled out due to the opposite motions 



of the ions in the electrolyte solution 



iii,y, 



10| . Therefore, in this study we neglected the hy- 



drodynamic interaction. Hydrodynamic interaction is important only when the chain length 



is very short H 



I2I. 

We varied the chain length (or the number of monomer) from 24 to 384 and studied 
the static and dynamic properties of PE under the action of an electric field, up to a field 
strength E = 2.0 ksT/lea). We set the monomer concentration Cm = O.OOOla"'^. In order 
to keep Cm constant, the size of the simulation box needs to change with A^. Instead of using 
a cubic simulation box, we chose a rectangular parallelepiped of 1.6A^cr x 79.06o" x 79.06cr, 
where the box size in the field direction is linearly proportional to to prevent overlap under 
periodic boundary condition when the chain unfolds. The added salt concentration was fixed 
at Cs = Cm/'i, the equivalence point C*. It has been shown that at this salt concentration. 
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FIG. 1: Re/Lc as a function of E at Cs = C* for different chain length N. The symbol 'x' denotes 
the inflection point of curve. 

the chain collapsed into a compact globule structure, in the absence of electric field Q, 
with its effective chain charge almost being neutralized. We performed Langevin dynamics 
simulations [isl with integrating time step equal to A t = O.OOSr. We ran firstly 10^ to 10^ 
time steps to bring the system to a steady state and then ran 10^ time steps to cumulate 
data for analysis. To simplify the notation, we assign in the following text that a, m, and 
ksT are the unit of length, mass, and energy, respectively. Therefore, the concentration will 
be described in unit of a~^, the strength of electric field in unit of ksT/^ea), and so forth. 



III. RESULTS AND DISCUSSIONS 



A. Degree of unfolding 

We start from studying the chain conformation under the action of an electric field. The 
degree of unfolding, defined as the ratio of the end-to-end distance Re of chain over the 
chain contour length Lc = {N — 1)6, is used to characterized the conformation. The results 
are plotted in Fig. [1] function of E. 
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Each curve in the plot denotes the variation of Re/Lc for a given chain length A^. We can 
see that when the electric field is weak, the ratio is a constant. This indicates an unperturbed 
conformation of chain and the chain remains in a collapsed structure. An abrupt increase 
appears when E is increased over some critical value E* . Rf. can become as large as 90% of 
Lc if the applied field is very strong. This indicates a structural transition from a collapsed 
structure to an elongated structure. We noticed that the value of E* depends on the chain 
length. The longer the chain length, the smaller the E* will be. Moreover, this structural 
transition happens in an interval of E. The size of the interval decreases with increasing chain 
length. Although the transition becomes sharper when chain length is long, Rf. increases in 
a continuous way with E, which suggests a second-order transition. 



B. Critical electric field E* 



The dependence of the critical electric field E* on chain length A^ has been investigated 
in salt-free js, 6] and in tetravalent salt solutions j?]. Both of these studies showed that E* 
scales roughly as N~^-^ to unfold a condensed chain. It is now important to know if this 
scaling law is valid for a PE chain condensed by trivalent counterions. To verify it, we follow 
firstly the method proposed by Netz 5, ^: E* is calculated by equating the polarization 
energy Upoi = p- E/2 and the thermal fluctuation energy k^T . Here p is the dipole moment 
of the PE-ion complex induced by the electric field and calculated by p = Zieifi — rem) 
where is the position vector, running over all the particles inside the complex, and Tcm 
is the center of mass of the PE. The complex is considered as a set of particles, including 
monomers and ions, inside the region of a worm-shaped tube which is the union of the 
jointed spheres of radius = 3, centered at each monomer center. The component of p at 
the field direction, p^., is plotted against the field strength E in Fig. [2J 

As seen in the log-log plot, px increases linearly with E with a slope equal to 1, when 
E is small. This is the well-known linear response of a dielectric object, p^ = aE, which 
has been reported in the previous studies 

a a. 

But different to the previous, we found 
that this linear region terminates before intersecting with the dotted line which denotes 
the relation pxE/2 = ksT , specially when the chain length is long. This is simply because 
the binding force to condense the PE chain in the trivalent salt solutions is weaker than 
in the tetravalent salt [7]. For the system studied by Netz , the chains were strongly 
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FIG. function of E for different chain length N. The dotted Une denotes the equation 

condensed because of the un-realistically strong Coulomb coupling chosen by him. Therefore, 
his method can be used only as a rough estimation of E* for the case of strong condensation 
but not suitable for the weak condensation. If we continue going with his method and 
calculate the intersection between the extended linear region and the dotted line, we will 
find that E* scales as A^-o-463{4) ^g^g open circles in Fig. [3t^a)). This scaling law seems to 
follow the prediction of Netz, N~^'^^'^, because the chain swelling exponent in zero electric 
field is z/ = 0.321(2) for this case (cf. Fig. [3t^b)). Nevertheless, E* obtained by this method 
is actually overestimated, going much over away the linear response region, specially when 
the chain is long. 

To give a more accurate estimation of E*, we follow here the definition of the unfolding 
electric field by taking simply the electric field at the inflection point of the curve Re/Lc 
vs. E. The inflection point on each curve in Fig. [1] is indicated by the symbol 'x'. The 
scahng law obtained by this method reads as ~ N^o.77ii) j^g^g close squares in Fig. [3](a)). 
The exponent —0.77(1) is significantly smaller than the one obtained by the Netz's method. 
Therefore, E* is smaller than Netz's prediction for a long chain and it is easier to unfold a 



7 




100 100 

N N 

FIG. 3: (a) E* vs. where the open circles denote the data obtained by Netz's method and the 
close squares denote the ones obtained from the inflection points, (b) Radius of gyration Rg vs. 
in zero electric field. 

PE chain in trivalent salt solutions. For example, for a chain of length 10^, E* is 1.76 x 10~^ 
according to the scaling law. This E* corresponds to about 185 V/cm, much smaller than 
2kV/cm predicted for the chains condensed by tetravalent salt in simulations Our results 
show that the valence of the condensing agent plays an important role in determination of 
the scaling law. There must exists more complicated mechanism to polarize and to unfold 
a PE chain in an electric field than our thinking. This mechanism will be investigated in 
detail in the future. 

C. Electrophoretic mobility and ion condensation 

We now study the electrophoretic mobility fipe of PE chain in electric fields of different 
strength and show how fipe changes with E when the chain is unfolded to an elongated 
structure, /ipe was calculated by Vpe/ E where Vpe is the velocity of the center of mass of the 
chain in the field direction. The results are shown in Fig. HI 

In weak electric fields, /ipe is nearly zero, indicating that the PE chain is effectively 
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charge- neutral, as reported in experiments [l^. While E is increased over E*, fipe turns to 
be negative and the chain starts to drift opposite the field direction, which suggests a negative 
effective chain charge. We found that the stronger the field, the faster the chain will drift. 
For a long chain, fip^ shows furthermore a plateau region when E > E*. The dependence 
of fipe on the electric field and the chain length gives a plausible way to electrophoretically 
separate PE chains by size in free solutions by means of chain unfolding transition [7|. 

The variation of be related to the ion condensation on the chain under the action 

of the electric field. Therefore, we studied here the number of the condensed trivalent ions 
on the chain by counting the ions inside the worm-shaped tube of radius = 3 around the 
chain. The results for = 384 in different strength of electric field are plotted in Fig. [5] 
against the monomer index l, rescaled from to 1, where l = denotes the first monomer 
heading toward the field direction and l = 1 denotes the last monomer of the other chain 
end. 

We saw that Nc{l) is flat when the applied field is small, E < 0.007, which shows 
an uniform distribution of the condensed trivalent ions along the chain. There is about 
0.33 trivalent counterions condensed on each monomer, which indicates the neutralization 
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FIG. 5: Nc{i) obtained in different field strength. The value of E is indicated near each curve. 

of the negatively-charged chain backbone by these condensed counterions. If we further 
increases the electric field, these condensed ions distribute non-uniformly on the chain where 
fewer ions condensed near the heading end {l = 0) than the tailing end {l = 1). When 
0.07 < E < 0.2, Nc{l) looks similar to an inclined line and the slope increases with E, 
resulting in a decrease of the total number of the condensed trivalent counterions on the 
chain. Therefore, \npe\ increases with E due to this partial detachment of the condensed 
ions by the electric field. At this moment, the PE-ion complex is polarized in a way that 
the condensed trivalent counterions are bound, basically immobile, on the chain. For the 
higher electric field, 0.2 < E < 1.0, Nc{l.) becomes a horizontal sigmoidal curve and the 
value in the middle chain region is independent of E. The appearance of this horizontal 
region reflects the fact that the condensed trivalent counterions are now gliding on the 
chain. These ions can be stripped off the chain by the strong electric field and the other ions 
in the bulk solution then condense onto it, establishing a steady state. The total number of 
the condensed trivalent counterions is approximately a constant in this electric field, which 
results in the plateau region of against E. For an even stronger electric field, such as 
E = 2.0, the baseline of the horizontal sigmoidal curve moves downward. The condensed 
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FIG. 6: (Color on line) Snapshots of unfolded PE chains in electric fields. The yellow, the white, 
the red, and the green spheres represent, respectively, the monomers, the monovalent counterions, 
the trivalent counterions, and the coins. The chain length, the electric field, the field direction, 
and the chain drifting direction are indicated in the figure. 

trivalent ions are stripped off the chain even more. The effective chain charge is thus more 
negative and |/ipe| increases. 

D. Conformation of an unfolded PE chain 

In our simulations, the PE chains were unfolded, for the most of the time, to an extended 
structure, similar to a straight line, aligned parallel to the field direction (see in Fig. [6](a)). 
Nonetheless, we observed sometimes that they were unfolded to a U-shaped structure in the 
electric fields. The open side of the U shape can point opposite or toward the chain drifting 
direction as shown in Fig. [HI panel (b) and (c), respectively. 

This U-shaped structure has been observed experimentally in electrophoresis of micro- 
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FIG. 7: Nc{t) for the two U-shaped chains in Fig. [6|^b) and (c), respectively. 
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171] and also been shown in simulations of the elastic uncharged/charged chains in 



stokes flows or in electric flelds 
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3, 



2l| . These studies showed that a combination of 



the elastic and the hydrodynamic effect results in the bending of a rigid chain into a horse- 
shoe shape, oriented perpendicular to the direction of motion 18|, |21| . If chains are charged 
and the driving force is an electric fleld, other effect, the electric polarization of the PE 
complex, will play a role, which favors parallel orientation to the electric fleld, and compete 



with the elasto-hydrodynamic effect 



191]. In our simulations, the PE chains are flexible and 



the hydrodynamic interaction is neglected. Therefore, different mechanism drives the chains 
to form U-shaped structures where the fleld-induced dipole moments on the two branches 
of a U-shaped chain establish an equilibrium. This phenomenon can be seen by plotting in 
Fig. [7] the distribution of the condensed trivalent counterions Nc{l) for the two U-shaped 
chains from Fig. [6](b) and (c), respectively. 

The symmetry of Nc{l) with respect to the middle point of the chain {l = 0.5) shows that 
an equilibrated polarization was established on the two branches of the U-chain in the electric 
fleld. The existence of two pointing directions of the open side of the U chain is a feature 
specially for the electric polarization. It is distinguishable to the elasto-hydrodynamic effect 
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where only the U-shaped structure with the open side opposite to the moving direction is 
produced. Moreover, we notice that the electrophoretic mobihty of an U-shaped chain is 
approximately equal to that of an elongated chain of half of the chain length. For example, 
Upe is —0.209(4) in Fig. Mjo), close to the mobility of the elongated chain of = 96, 
—0.225(3). Furthermore, we verified the stability of these U-shaped chains and found that 
they can persist through the whole simulation period corresponding to, at least, the order 
of microsecond. However, by introducing some perturbations such as AC electric fields, the 
U-shaped structure can be transfered into the elongated chain structure but the inverted 
direction of transfer cannot be realized. Therefore, the U-shaped structure is probably 
metastable. We have calculated the total energy of the system for the U-shaped chain 
structure and also for the extended-chain structure. We found that the previous energy is, at 
least, 5% higher than the latter. Moreover, the U-shaped chain has a slower electrophoretic 
mobility than the extended chain, which implies a larger number of counterions condensed 
on the U-shaped chain to decrease the effective chain charge; consequently, fewer ions are 
presented in the bulk solution and the entropy of the solution is small, compared to the 
extended-chain structure. Therefore, the free energy of the system is lower for the extended 
chain than for the U-shaped chain. This estimation supports that the U-shaped structure 
is metastable. Since the open side of the U-shaped chain can point to one of the two 
directions, along or against the field direction, we predict the existence of other metastable 
states, due to polarization, in which the chain shows many bends, such as S-shaped or 
W-shaped structures, in electric fields. 



IV. CONCLUSIONS 



We have studied the behavior of single polyelectrolytes condensed by trivalent salt under 
the action of an uniform electric field by means of Langevin dynamics simulations. We found 
that the chains unfolded while the strength of the electric field is stronger than some critical 
value E*, similar to the previous study where the chains were condensed by tetravalent 
salt 0]. E* shows scaling- law dependence on the chain length A^, reading as E* ~ ]\[~^-'^'^(^) , 
The exponent in the scaling law is different from the prediction by Netz j^, 6] and from 
the simulations in tetravalent salt solutions j?], which demonstrated the importance of the 
salt valence on the exponent. Therefore, the weaker the condensing agent, the larger the 
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absolute value of the exponent and the easier the unfolding of a condensed chain will be. We 
showed that the electrophoretic mobility of chain \^ipe\ drastically increases while the chain 
is unfolded. The distribution of the condensed counterions on the chain was studied and 
related to the change of the mobility in different regions of electric field. The dependence 
of /ipe on the chain length and the electric field enables us to device a way to impart chain- 
length dependence in free-solution electrophoresis through chain-unfolding mechanism in 
electric fields. Finally, we pointed out the possibihty to unfold a condensed PE chain into 
an U-shaped structure in electric fields, in addition to the elongated structure, with the open 
side of the U heading or taihng the chain drifting direction. This structure is a result of 
purely electric polarization, different from the formation of the horseshoe-shaped chains in 
sedimentation experiments caused by elasto-hydrodynamics. 
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